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The  interference  modulator  (IM)  is  a  microwave  source  concept  that  can 
produce  a  broadband,  amplitude  modulated  radio  frequency  (rf)  waveform 
with  independently  tunable  carrier  and  modulation  frequencies.  In  its  simplest 
form,  the  IM  is  a  two-channel  device  that  can  produce  an  ampi  itude  modulated 
signal  at  50-percent  duty  by  either  generating  and  beating  two  electromagnetic 
waves  together  in  the  far  field  or  by  combining  the  two  signals  with  a  four-port 
hybrid  network.  With  additional  channels,  an  IM  can  be  made  to  synthesize 
more  complex  waveforms  and  duty  factors.  A  general  methodology  for  the 
design  of  an  n-channel  (n  >  2)  IM  is  developed.  As  a  validation  of  the  design 
methodology,  a  two-channel  prototype  IM  was  constructed.  The  design, 
fabrication,  and  testing  of  the  two-channel  prototype  are  described  in  detail. 
Experiments  carried  out  with  the  prototype  IM  have  successfully  demonstrated 
the  IM  concept,  producing  a  source  with  a  modulation  envelope  that  is 
continuouslytunableoverfourdecades  (from  100  Hz  to  1  MHz),  with  a  carrier 
frequency  that  is  continuously  tunable  from  800  MHz  to  2  GHz. 
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Introduction 


The  interference  modulator  (IM)  is  a  microwave  source  that  can  gener¬ 
ate  a  broadband,  amplitude-modulated  radio  frequency  (rf)  waveform 
with  independently  tunable  carrier  and  modulation  frequencies.  Two 
key  features  of  the  IM  are  (1)  cw  operation  of  the  rf  source  components, 
which  minimizes  the  power  conditioning  requirements,  and  (2)  opera¬ 
tion  over  a  very  broad  bandwidth.  These  features  make  the  IM  a  highly 
flexible  microwave  source  with  applications  in  studies  of  electromag¬ 
netic  (EM)  coupling  and  the  susceptibility  of  electronic  components 
and  systems.  This  report  presents  the  methodology  for  the  design  of 
an  IM  and  describes  experiments  carried  out  with  a  prototype. 

1.1  Motivation  for  Development  of  IM 

EM  radiation  may  couple  into  a  complex  object  through  a  variety  of 
paths.  Front-door  coupling  takes  place  when  an  externally  applied  sig¬ 
nal  enters  the  system  along  an  intentional  path  for  EM  radiation,  such 
as  an  antenna  or  sensor  aperture.  Back-door  couphng  takes  place  when 
the  radiation  enters  the  system  along  an  unintentional  path,  such  as 
cables  or  cracks  in  the  system  housing. 

Amplitude-modulated  rf  signals  can  be  more  eflFective  at  jamming 
or  upsetting  electronic  systems  than  non-amplitude-modulated  signals. 
Amplitude  modulation  could  contribute  to  effective  electronic  upset  in 
the  following  manner:  the  high-frequency  (short  wavelength)  carrier 
could  allow  the  rf  signal  to  enter  the  system  via  a  front-  or  back-door 
path,  such  as  a  sensor  port  or  a  seam  in  the  hardware  enclosure;  once 
inside  the  target  system,  the  lower  frequency  modulation  (if  properly 
selected)  could  be  rectified  by  nonlinear  electronic  circuit  elements,  such 
as  transistor  junctions,  and  be  erroneously  interpreted  as  valid  data  by 
the  system  logic.  If  this  occurs,  a  failure  or  upset  of  the  target  system 
may  result. 

Significant  changes  in  the  efficiency  of  rf  coupling  into  enclosures  can 
be  generated  by  small  changes  in  rf  frequency  or  modulation.  Coupling 
effectiveness  may  change  by  as  much  as  20  to  30  dB  for  changes  in 
frequency  as  small  as  a  few  megahertz.  Because  of  this  sensitivity, 
the  ideal  source  for  evaluating  the  rf  susceptibUity  of  a  given  system 
and/or  potentially  upsetting  it  should  have  a  carrier  frequency  and 
modulation  envelope  with  independent  and  near-infinite  tunabilities. 
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Commercially  available  components  do  not  have  the  desired  bandwidth 
and/or  tunability  for  these  applications.  The  IM  is  an  attempt  at  a 
practical  design  of  a  device  with  a  widely  selectable  carrier  frequency 
and  independently  tunable  modulation  envelope. 

1.2  IM  Basics 

Two  examples  of  IM  configurations  are  shown  conceptually  in  figure  1. 
A  high-frequency  carrier  signal,  /c,  and  a  lower  firequency  modulation 
signal,  fm,  are  the  inputs  to  the  IM.  For  outputs,  the  IM  will  pass  the 
carrier  frequency,  fc,  and  generate  a  processed  signal  consisting  of  the 
carrier  frequency  plus  (or  minus)  the  modulation  frequency,  fc+fm  {fc- 
fm)-  These  two  output  signals  can  be  post-amplified  and  radiated  into 
free  space  through  a  pair  of  matched  antennas.  In  the  far  field,  fc  and 
fc  +  fm  {fc  —  fm)  will  linearly  add  to  produce  an  amplitude-modulated 
EM  wave  with  a  50-percent  duty  factor  and  a  peak  ampUtude  that  is 
twice  that  of  the  individual  waves. 

To  illustrate  the  modulation  effect  of  linearly  combining  two  waves 
closely  spaced  in  frequency,  figure  2  presents  a  plot  of  two  sinusoidal 
waves  with  equal  amplitudes,  one  with  a  frequency  of  1.0  GHz,  and 
the  other  with  a  frequency  of  1.1  GHz.  The  sum  of  the  two  waveforms 
is  plotted  in  figmre  3,  demonstrating  the  envelope  modulation  of  the 
carrier  frequency.  Note  that  the  peak  amplitude  of  the  sum  is  twice  that 
of  the  constituent  waveforms.  The  spectral  content  of  this  “beat  wave” 
is  a  superposition  of  the  frequency  spectra  of  each  of  the  constituent 
waves,  as  shown  in  figure  4,  the  Fourier  transform  of  figure  3.  The 
broadening  of  the  spectrum  in  figure  4  is  a  consequence  of  the  finite 
time  duration  of  the  wave  used  in  the  numerical  analysis. 

There  are  several  advantages  to  producing  an  amplitude-modulated 
signal  in  this  fashion,  including  minimization  of  the  system  volmne  and 
weight,  reduced  stress  on  individual  system  components,  and  modular- 
By  generating  the  modulation  envelope  through  the  interference 
of  two  waves  in  the  far  field  (see  fig.  1(b)  and  (c)),  we  can  operate 
the  two  rf  sources  in  a  cw  mode  where  the  individual  sources  can  op¬ 
erate  at  one-half  the  average  power  (one-fourth  the  peak  power)  of 
the  equivalent  single-antenna,  amplitude-modulated  source.  The  lower 
power  requirements  of  the  IM  concept  result  in  reduced  voltage  stress 
on  the  individual  components  (amplifiers,  antenna  interfaces,  transmis- 
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Figure  1:  Conventional  and  IM-based  schemes  for  producing  a  50- 
percent  duty,  amplitude-modulated  signal  in  far  field:  (a)  conventional 
single  antenna  and  pulsed  amplifier  {Ae  =  ^o);  (b)  two-channel  IM  sys¬ 
tem  (Ae  =  Ao);  (c)  two-element  phased-array  IM  system  {A^  =  Aq/2). 
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Frequency  (GHz) 


Figure  4;  Fourier  transform  of  beat  wave  of  figure  3. 

sion  systems)  as  well  as  a  reduced  volume  and  weight  compared  to  the 
equivalent  single  pulsed  rf  source  (depicted  in  fig.  1(a)).  Note  that  the 
equivalent  antenna  aperture,  Ag,  for  each  scenario  in  figure  1  has  been 
normalized  to  the  single  antenna  case  (fig.  1(a)).  Using  the  IM  concept, 
we  can  produce  the  far-field  wave  interference  effect  either  by  a  two- 
channel,  non-phase-synchronized  system  requiring  two  antennas  each 
of  aperture  .do  as  shown  in  figure  1(b),  or  by  a  two-element  phased 
array  requiring  two  antennas  each  with  aperture  Ao/2,  as  shown  in 
figure  1(c). 

In  addition,  the  IM  concept  can  be  extended  to  n  sources  (n  >  2); 
with  proper  phase  control  of  the  individual  sources,  more  complicated 
rf  waveforms  and  lower  duty  factors  can  be  synthesized.  Like  the  two- 
source  IM,  multiple-source  IM  systems  have  the  advantages  of  lower 
voltage  stress  on  the  individual  components  and  potentially  lower  sys¬ 
tem  weight  and  volume. 
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Mixer  1 


2  Development  of  a  Prototype  Interference  Modulator 

A  prototype  IM  was  developed  from  a  custom-designed  active-network 
single-sideband  modulator  and  commercially  available  microwave  com¬ 
ponents  (mixers,  hybrids,  and  power  combiner  and  dividers).  A  block 
diagram  of  the  IM  is  shown  in  figme  5.  The  IM  carrier  frequency,  /c, 
is  synthesized  from  two  components  —  a  constant  upshift  frequency, 
/„/,  and  a  tunable  component,  fr/,  such  that  fc  =  fuf  +  frf-  The  IM 
modulation  frequency,  fcm,  is  synthesized  from  frf  and  a  tunable  low- 
frequency  component,  fm,  which  has  been  frequency  upshifted  by  /„/ 
so  that  fcm  =  (frf  +  fuf)  +  fm  =  fc  +  fm-  For  the  highest  generation 
efficiency,  /„/  -I-  fm  is  produced  by  a  single-sideband  network  composed 
of  a  broadband  90°  phase-difference  network  with  a  power  combiner 
and  two  double-balanced  mixers  to  perform  the  frequency  upshift. 

2.1  Single-Sideband  Generation  Using  1 8z  Q  Modulated  Signals 

2.1.1  Basic  Approach 

An  /  &  (^  modulation  network  takes  a  single  sinusoidal  waveform  as 
input  and  generates  two  signals  as  output:  an  in-phase  signal  (7)  and 
a  quadrature  signal  ( Q)  that  is  90°  out  of  phase  with  the  input  signal. 
Such  a  network  is  also  referred  to  as  a  90°  phase-difference  network. 
The  I  and  Q  outputs  can  be  further  combined  with  microwave  mixers, 
hybrids,  and  combiners  to  generate  a  single-sideband  signal. 
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Figure  6:  Block  diagram  of  IM  single-sideband  modulation  subsystem. 


Figure  6  is  a  block  diagram  of  the  single-sideband  modulation  sub¬ 
network  of  the  IM.  The  input  waveform  to  the  90°  phase-difference 
network  has  the  form 


fm{0J,t)  =  Vm  sin{Umt),  (1) 

where  Vm  is  the  voltage  amplitude  in  volts,  Um  is  the  modulation  fre¬ 
quency  in  radians  per  second,  and  t  is  the  time  in  seconds. 

After  processing  by  the  90°  phase-difference  network,  the  I  and  Q 
output  signals  have  the  form 


=  Vml  Sin{ujmt), 

(2) 

t)  =  VmQ  sin{Umt  -H  (f>m), 

(3) 

where  Vmi  and  VmQ  are  the  voltage  amplitudes  of  the  I  and  Q  signals, 
respectively,  and  (f>m  is  the  phase  difference  introduced  by  the  network. 
Ideally,  4>m  =  90°,  but  we  treat  it  as  a  variable  in  order  to  perform  an 
error  sensitivity  analysis. 

The  /  and  Q  signals  are  then  mixed  with  a  sinusoidal  upshift  signal, 
/„/,  as  shown  in  figure  6.  The  input  upshift  signal  is  described  by 

fuf{u},t)  =  Vuf  siniuJuft),  (4) 

where  K/  is  the  voltage  amplitude  of  the  upshift  signal  and  ujuf  is  the 
upshift  frequency  in  radians  per  second. 
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After  passing  through  a  0°  power  divider  and  a  3-dB  90°  hybrid, 
the  upshift  signal  is  split  into  in-phase  and  quadratiure  signals: 

fuf,i  =  ^  sin(a;„/t),  (5) 

fuf,Q  =  ^  siniujuft  +  |).  (6) 

As  shown  in  figure  6,  the  waveforms  of  equations  (2)  and  (5)  are 
combined  in  Mixer  1,  and  the  waveforms  of  equations  (3)  and  (6)  are 
combined  in  Mixer  2.  Using  the  trigonometric  identities, 

sin  A  sin  =  -  [cos(A  -  B)  -  cos(A  -I-  B)] ,  (7) 

cos(A±^)  =  =FsinA,  (8) 

we  obtain  the  rf  outputs  of  Mixers  1  and  2,  respectively. 


Vmi  sin(a;,n^)  (^)  sin{ujuft  +  (9) 


1  V„,V.t  \ 

^cos 

(iOm  2 

1 

i 

5  v„,v„,  \ 

|cos 

[CJm  +  +  2 

) 

g  VmiVuf  [sin(a;  m  uJuf)t  +  sin(t(;r„  -f-  u)uf)t]  ■ 


/mix2(^)t)  —  KnQ  sin(a;m^  +  0m)  sin(a;u^t -f  — )  (10) 

g  ^mQ^uf  {cos  [(u^m  ^uf)t  0m] }  “I" 

g  ^mQ^uf  {cos  [{oJm  +  -|-  0m]}  • 

If  (f)m  —  7r/2  in  equation  (10),  then  equation  (10)  reduces  to 

/mix2(0m  —  '^/2)  =  —  [—  sin(a;,7i  —  Uuf)t]  -(- 

1 


(11) 


Combining  equations  (9)  and  (11)  using  the  0°  power  combiner 
shown  in  figure  6,  we  obtain 


/sum(t^,i)  =  ^  {{VraI-VmQ)sm{u;m-(^uf)t]+  (12) 

^  [{Vmi  +  VmQ)  sin(a;m  +  a;„/)t]  • 

If  Vmi  =  VmQ  =  Kio,  we  obtain  the  upper  sideband  signal,  /usb. 


/usb(t^,  ^)  =  \  ^moVuf  Sin(a;^  +  0Juf)t.  (13) 

Note  that  if  the  LO  inputs  to  Mixers  1  and  2  are  exchanged  (fig.  6), 
the  output  of  the  0°  power  combiner  will  produce  the  lower  sideband 
signal,  /isb'  ^ 

^)  =  4  VmoVuf  sin(a;„i  -  a;„/)i.  (14) 


2.1.2  Effect  of  Amplitude  and  Phase  Error  on  Single-Sideband  Gen¬ 
eration 

If  the  I  and  Q  signals  produced  by  the  90°  phase-difference  network  are 
not  of  equal  amplitudes  or  if  the  phase  difference  deviates  from  90°, 
then  when  the  I  and  Q  signals  are  summed  in  the  0°  power  combiner 
(fig.  6),  the  unwanted  (upper  or  lower)  sideband  will  not  be  completely 
suppressed.  To  evaluate  the  effect  of  amplitude  and  phase  error  on 
single-sideband  generation,  we  first  sum  equations  (9)  and  (10): 

/mixi  +  /mix2  =  ^  K/  {[Vmi  “  VmQ  sixi (f)m]  sin(a;„,  -  u;„/)t  (15) 

+  [Vmi  +  VmQ  sin  4)m]  sin(a;,„  +  a;„/)t 
-f  [1  -  COs(a;r„  +  <^uf)t]  cos  <l)m  COs{uJm  “  Wu/)0  • 


If  the  /  and  Q  phase  difference,  (f>mi  is  close  to  90°,  cos  (pm  —  0 
and  the  last  line  of  equation  (15)  can  be  neglected.  The  remaining 
two  lines  of  equation  (15)  represent  the  lower  (LSB)  and  upper  (USB) 
sideband  terms,  respectively.  Taking  the  ratio  of  their  amplitudes  yields 
the  relative  suppression  (enhancement)  of  one  sideband  over  the  other. 
From  equation  (15), 


LSB 

USB 


=  20  log 


Vmi  -  VmQ  sin  <P 

m 

Vmi  +  VmQ  sin  <pm 


(16) 
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R4 


As  an  example,  let  VmQ  =  0.814»/  and  =  88°.  Then  the  suppres¬ 
sion  of  the  lower  sideband  in  decibels  would  be 

T  SR 

USB  =  “ 

Thus,  a  20-percent  diflPerence  in  the  I  and  Q  voltage  ampUtudes  and  a 
<3-percent  phase  error  combine  to  produce  a  nonideal  sideband  sup¬ 
pression  of  —19.06  dB. 

2.2  Synthesis  of  a  Broadband  90°  Phase-Difference  Network 

One  of  the  unique  components  of  the  IM  design  is  a  broadband  90° 
phase-difference  network  capable  of  operating  over  a  frequency  range 
of  four  decades.  The  basic  building  block  of  the  circuit  is  an  active  all¬ 
pass  filter  network,  depicted  in  figure  7.  Using  the  synthesis  method 
described  by  Sebol,^  we  can  use  all-pass  filter  blocks  to  realize  an  even- 
order,  broadband,  active  network  that  takes  a  single  frequency  as  input 
and  produces  in-phase  (I)  and  quadratme  (Q,  i.e.,  90°  out  of  phase) 
signals  as  outputs. 

The  synthesis  of  the  90°  phase-difference  network  begins  with  the 
specification  of  the  frequency  range  over  which  the  network  is  to  operate 

^Ronald  Sebol,  “Design  of  Active  90-Degree  Phase  Difference  Networks,”  Harry 
Diamond  Laboratories,  HDL-TM-68-18  (May  1968). 
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and  a  maximum  phase  error  tolerance,  6,  in  radians.  The  frequency 
range  is  bounded  by  a  starting  frequency,  and  a  cutoff  frequency, 
fb-  Prom  Bedrosian,^  the  design  equations  for  the  filter  network  are 
given  by 


U  _ 

(18) 

II 

1 

(19) 

K  =  Kik), 

(20) 

K'  =  K{k'), 

(21) 

1 

11 

(22) 

where  K  is  the  complete  elliptic  integral  of  the  first  kind  and  n  is  the 
order  of  the  filter. 

Prom  the  above  equations,  we  can  determine  the  required  number  of 
filter  sections  (poles  or  zeros),  n,  for  a  given  frequency  band  and  phase 
ripple.  Once  the  munber  of  poles  (zeros)  is  determined,  the  location  of 
the  pole  (zero)  is  calculated  from  the  elliptic  functions  cn  and  sn 
with  the  formula 


\fa)  sn{uj,k) 

where 

J=0.1.2,....n-1.  (24) 

Using  the  relationships  between  the  elliptic  functions  and  the  theta 
function,  we  can  express  equation  (23)  in  a  series  expansion^ 

^S.  D.  Bedrosian,  “Normalized  Design  of  90®  Phase-Difference  Networks,”  IRE 
Trans.  Circuit  Theory  (June  1960),  128-136. 

^Milton  Abramowitz  and  Irene  Stegun,  eds..  Handbook  of  Mathematical  Func¬ 
tions^  New  York:  Dover  Publications,  Inc.,  1972,  p.  575. 
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(25) 


y)  cos  (21  +  1)  Tj 

„  _  ^=0 _ 

Pj  oo  ? 

£(_l)'g'a+i)siii(2Z  +  l)r,- 

/=o 

where 

r,  =  (-l)^(2j  +  l)^,  j=0,l,2,...,n-l,  (26) 

q  =  e~^^.  (27) 

A  C++  program  was  written  to  evaluate  equation  (25)  and  to  cal¬ 
culate  the  poles  (zeros)  and  response  (phase  error)  of  the  filter  as  a 
function  of  /„,  ft,  and  n.  Appendix  A  contains  a  listing  of  this  pro¬ 
gram. 

Once  the  network  order  has  been  determined  by  the  phase  error 
tolerance  (maximum  allowable  deviation  from  90°)  and  the  bandwidth 
ratio,  fb/fa,  and  the  pole  (zero)  pairs  have  been  generated  by  the  C++ 
program,  the  component  values  for  each  stage  of  the  filter  network  can 
be  calculated.  Referring  to  the  nomenclature  of  figure  7  and  using 
Sebol’s  synthesis  equations,^  we  compute  the  component  values  as 


P  (Pa+  Pb)  ±  QPaPb  +  Pb 

2paPb 

(28) 

PaPbR2  ’ 

(29) 

..  2^2 

nr 

(30) 

i?4  =  KR3, 

(31) 

where  Pa  and  pb  are  pole  pairs  computed  from  equation  (25).  If  we 
choose  i?3  =  1  kfl,  then  equations  (28)  to  (31)  are  sufficient  to  synthe¬ 
size  the  filter.  Note  that  the  poles  must  be  paired  so  that  the  discrimi¬ 
nant  of  equation  (28)  is  always  greater  than  or  equal  to  zero.  Appendix 
B  contains  a  listing  of  a  C  program  that  will  automatically  generate 
the  network  component  values. 
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Table  1:  Design  specifications  for  90°  phase-difference  network. 


Specification 

Value 

Frequency  range 

100  Hz  <  /„!  <  1  MHz 

Phase-error  tolerance 

±1° 

Carrier/sideband  rejection 

>  25  dB 

3  Design  Example 

As  a  design  exercise,  we  developed  an  IM  with  a  carrier  frequency 
that  was  continuously  tunable  over  a  frequency  range  of  approximately 
800  MHz  to  2  GHz  and  a  modulation  frequency  continuously  tunable 
from  100  Hz  to  1  MHz.  Table  1  summarizes  the  specifications  for  the 
90°  phase-difference  network. 

Prom  table  1,  the  bandwidth  ratio,  defined  as  the  ratio  of  the  upper 
to  the  lower  frequency,  is  given  by  fb/ fa  =  10, 000.  Using  equation  (22) 
and  the  code  PHSESHFT.CC  (app  A),  we  found  the  network  order  to 
be  n  =  12.  Assuming  a  phase  error  tolerance  of  ±1°,  we  calculated  the 
roots  from  equation  (25)  (using  the  code  PHSESHFT.CC). 

The  optimum  pole  pair  combinations  for  the  active  network  were 
determined  with  the  code  IQ_SYNTH.C  (app  B),  which  also  generated 
the  component  values  for  each  operational  amplifier  stage. 

3.1  Design  £uad  Fabrication  of  90°  Phase-Diflference  Network 

As  synthesized  by  the  method  described  in  section  2.2,  the  90°  phase- 
difference  network  used  six  high-gain/bandwidth-product  operational 
amplifiers  to  realize  the  in-phase  and  quadrature  outputs  over  the  re¬ 
quired  bandwidth.  Several  SPICE  circuit  models  of  the  90°  phase- 
difference  network  were  developed  with  both  ideal  and  manufactiurer- 
supplied  operational  amplifier  models.  We  modified  the  ideal  (com¬ 
puted)  component  values  in  the  SPICE  models  to  correspond  with  the 
values  of  commercially  available  components  to  ensure  that  the  effect 
of  the  nonideal  components  on  circuit  performance  was  negligible.  The 
schematic  for  the  network  is  shown  in  figure  8.  Based  on  the  results  of 
the  SPICE  modeling,  a  printed  circuit  board  was  fabricated,  populated 
with  components,  debugged,  and  tested. 
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R14  R24  R34 


Figiire  8;  90°  phase  difference  network  schematic. 


3.1.1  Printed  Circuit  Board  Layout  and  Fabrication 

The  mask  for  a  two-layer  printed  circuit  board  (PCB)  was  generated 
on  a  DEC  workstation  with  the  aid  of  two  CAD/CAM  programs.  Pow- 
erview  (Viewlogic)  was  used  to  generate  a  circuit  schematic  and  nethst 
from  the  SPICE  input  file.  The  netlist  contained  a  logical  map  of  the 
circuit  components  and  their  associated  pin  connections.  SCICARDS 
(Harris)  used  this  information  to  generate  the  physical  layout  of  the 
PCB,  including  all  the  component  locations,  trace  routings,  and  vias. 

The  two-layer  PCB  measured  approximately  15  x  23  cm  (6  x  9  in.). 
The  board  design  includes  electrical  test  points  at  each  operational- 
amplifier  stage  to  facilitate  debugging.  Board  power  can  be  provided 
externally  from  a  dc  power  supply  or  internally  from  a  set  of  four  9-V 
batteries.  Gold-plated,  right-angle  SMA  connectors  were  moimted  on 
the  board  to  accept  the  modulation  frequency  as  input  and  to  provide 
access  to  the  in-phase  (I)  and  quadrature  (Q)  outputs  of  the  network. 
A  tunable  resistive  divider  network  consisting  of  a  series  10-turn  200- 
Q,  potentiometer  and  a  shunt  51-0  resistor  was  placed  at  both  the  I 
and  the  Q  output  nodes  so  that  their  respective  amplitudes  could  be 
matched. 
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All  the  components  used  standard  library  shapes,  except  for  the 
electrolytic  bypass  capacitors,  the  battery  holders,  and  the  right-angle 
SMA  receptacles.  These  shapes  were  defined  in  the  user  hbrary.  The 
hole-to-pad  ratio  for  the  user-defined  shapes  was  approximately  two 
to  one.  This  ratio  provided  a  sufiicient  amount  of  contact  area  when 
the  discrete  components  were  soldered  onto  the  board.  Surface-mount 
resistors  (1-percent  tolerance,  size  1206)  and  capacitors  (10-percent  tol¬ 
erance,  size  1206)  were  used  except  for  the  bypass  capacitors.  SPICE 
simulations  indicated  that  the  circuit  performance  was  highly  sensitive 
to  the  capacitor  values:  a  10-percent  deviation  in  the  ideal  capacitance 
value  could  produce  undesirable  results.  Therefore,  we  hand-selected 
each  chip  capacitor  to  ensure  that  it  was  as  close  as  possible  to  its  ideal 
value;  in  addition,  the  circuit  board  contained  extra  pads  in  parallel 
with  each  capacitor  to  allow  for  additional  trimming  if  necessary. 

Although  SCICARDS  automatically  generated  a  physical  layout 
from  the  netlist,  we  edited  this  layout  extensively  to  ensure  optimiun 
trace  widths,  routing,  and  positioning  of  all  components.  The  criteria 
used  for  optimizing  the  trace  routing  were  to  minimize  electrical  path 
lengths  and  vias,  and  to  simplify  the  board  layout  so  as  to  facilitate 
assembly  and  testing.  The  traces  were  kept  a  reasonable  distance  firom 
any  holes  or  solder  points  on  the  board  so  that  components  would  not 
be  soldered  to  an  incorrect  path  or  pin.  Line  segments  were  reduced 
and  kept  as  straight  and  long  as  possible.  Unnecessary  vias  were  elimi¬ 
nated,  and  vias  placed  too  close  to  components  or  nets  were  relocated. 
The  default  via  width  was  set  at  0.062  in.  We  widened  the  power  sup¬ 
ply  and  ground  trace  widths  with  respect  to  the  other  traces  to  reduce 
the  resistance  in  the  circuit  (0.004  in.  as  compared  to  0.002  in.  for  non¬ 
supply  or  ground  circuit  traces).  Test  points,  designed  as  a  small  hole 
within  a  large  pad,  were  included  at  the  input  and  output  of  the  each 
amplifier  stage. 

The  board  used  six  very-high-speed  operational  amplifiers  (Analog 
Devices  AD846,  small  signal  bandwidth  =  46  MHz).  Each  operational 
amplifier  used  two  bypass  capacitors  (a  2,2-/iF-  electrolytic  and  a  ce¬ 
ramic  0.1  /xF)  in  parallel  across  the  plus  and  minus  voltage  terminals 
to  filter  out  unwanted  frequencies  generated  by  the  power  supply.  If  left 
unfiltered,  these  frequencies  could  lead  to  instabilities  in  the  operation 
of  the  circuit.  The  bypass  capacitors  were  placed  as  close  as  possible 
to  the  operational  amplifiers  to  minimize  stray  inductances. 
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Figure  9:  Typical  oscilloscope  trace  of  /  and  Q  outputs,  demonstrating 
90°  phase  shift. 

Once  fully  populated,  the  PCB  was  mounted  in  a  standard  19-in. 
rack  mount  chassis  along  with  the  microwave  mixers,  hybrids,  and  com¬ 
biners  that  completed  the  IM  system.  A  47- /jF  electrol3d;ic  capacitor 
and  a  0.1-/iF  ceramic  capacitor  were  used  to  bypass  the  external  power 
supply,  and  a  0.4-/zF  capacitor  was  used  to  bypass  the  battery  power 
supply  mounted  on  the  board. 

3.1.2  90°  Phase-Difference  Circuit:  Measured  Performance 

When  first  fabricated,  the  90°  phase-difference  board  exhibited  high- 
frequency  oscillations  in  both  the  /and  Q  output  signals.  We  eliminated 
this  problem  by  placing  two  small- valued  capacitors  across  the  inverting 
and  non-inverting  inputs  of  the  third  stage  (I)  op-amp  and  by  placing 
one  small- valued  capacitor  across  the  inverting  and  noninverting  inputs 
of  the  sixth  stage  (Q)  op-amp.  We  matched  the  output  amplitudes 
of  the  I  and  Q  signals  using  the  potentiometers  located  at  the  two 
outputs.  The  output  amplitudes  of  both  branches  remained  matched 
and  constant  over  the  designed  frequency  range  of  100  Hz  to  1  MHz. 
Figure  9  is  a  typical  time-domain  oscilloscope  trace  of  the  /  and  Q 
outputs  demonstrating  the  90°  phase  shift. 

The  phase  difference  between  the  /  and  Q  outputs  as  a  function 
of  frequency  was  measured  with  an  HP  3575A  gain-phase  meter.  The 
experimentally  measured  phase  difference  as  a  function  of  modulation 
frequency  is  plotted  in  figure  10,  along  with  the  SPICE  predictions 
using  both  an  ideal  op-amp  model  and  an  Analog  Devices  AD846  op- 
amp  model. 

The  corresponding  measured  and  calculated  phase  error  (i.e.,  devi- 


Figure  10:  Calculated  and  measured  phase  difference  between  I  and  Q 
outputs  of  90°  phase-difference  network. 

ation  from  90°)  is  plotted  in  figure  11.  Using  ideal  components  (op- 
amps  with  infinite  gain-bandwidth  products),  the  SPICE  simulations 
show  that  the  circuit  synthesis  method  had  achieved  its  design  goals  of 
±1°  phase  error  tolerance  over  a  frequency  range  of  100  Hz  to  1  MHz. 
The  predicted  performance  at  the  lower  frequencies  (<10  kHz)  suffers 
when  the  more  realistic,  manufactured-supplied  AD846  op-amp  model 
is  used.  The  measiured  phase  error  is  poorer  than  either  the  ideal  or 
nonideal  simulations,  probably  because  of  stray  capacitances  and  induc¬ 
tances  that  were  not  taken  into  account  in  the  SPICE  model.  However, 
as  discussed  in  section  3.2,  the  experimentally  achieved  phase  error  was 
more  than  sufficient  for  the  performance  of  the  prototype  IM. 

3.2  Single-Sidebeind  Modulator  Performance 

After  debugging,  the  90°  phase-difference  network  board  was  integrated 
into  the  single-sideband  (SSB)  modulator  subsystem  of  the  IM,  as 
shown  in  the  block  diagram  of  figure  6.  The  modulating  frequency, 
fm,  was  provided  by  an  HP651A  oscillator,  and  the  upshift  frequency, 
/„/,  was  provided  by  an  HP8656A  signal  generator.  Figure  12  is  a 
typical  plot  of  the  spectral  output  of  the  SSB  subsystem.  As  can  be 
seen  in  the  figure,  the  carrier  and  unwanted  sideband  frequencies  are 
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Figure  11:  Calculated  and  measured  phase  error  of  the  90°  phase- 
difference  network. 


suppressed  at  least  —25  dB. 

A  key  requirement  for  the  SSB  subsystem  was  operation  over  a 
very  broad  range  of  frequencies.  The  broadband  behavior  of  the  SSB 
subsystem  is  plotted  in  figure  13,  where  the  measured  suppression  of  the 
carrier  and  imwanted  sideband  is  recorded  as  a  function  of  modulating 
frequency. 

As  evidenced  by  the  measured  data,  the  single-sideband  modulator 
design,  as  synthesized  by  the  techniques  of  sections  2.1  and  2.2,  is 
highly  successful,  producing  a  relatively  pure  tone  with  at  least  -25  dB 
of  suppression  of  unwanted  frequencies  over  a  frequency  span  of  four 
decades.  Over  the  design  frequency  range  of  100  Hz  to  1  MHz,  the 
measured  average  carrier  suppression  was  —25.5  dB,  and  the  measured 
unwanted  sideband  suppression  ranged  from  —26  to  —46  dB. 

3.3  Measured  Far-Field  Performance  of  Interference  Modulator 

To  evaluate  the  performance  of  the  IM  system,  we  post-amplified  the  IM 
carrier  (/c)  and  modulation  (/c  —  fm)  outputs  and  used  them  to  drive 
two  matched  standard-gain  horn  antennas.  Figure  14  is  a  schematic 
of  the  test  setup.  The  initial  outputs  from  the  IM  were  amplified  by 
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Figure  12:  Typical  measured  spectrum  of  single-sideband  modulator 
subsystem. 
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Figure  13:  Carrier  and  unwanted  sideband  suppression. 
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Figure  14:  Block  diagram  of  IM  system  test. 


Receive  hom 


a  pair  of  broadband  solid-state  amplifiers  (Mini-Circuits  ZFL-llAD, 
8-dB  gain).  The  outputs  from  the  solid  state  amplifiers  were  then  fed 
into  a  pair  of  Varian  20- W  traveling- wave  tube  amplifiers,  which  in  turn 
were  connected  to  a  pair  of  matched  WR650  waveguide  horns  spaced 
43.2  cm  apart  (edge-to-edge  spacing).  The  horns  radiated  into  a  large 
anechoic  chamber.  A  third  WR650  horn  was  placed  in  the  far  field  at 
a  distance  of  4.3  m  from  the  transmitting  horns  to  measure  the  output 
radiation. 

Initially,  the  time-  and  frequency-domain  output  behavior  of  each 
of  the  transmitting  horns  was  evaluated  separately.  With  the  IM  set  up 
to  suppress  the  upper  sideband,  far-field  measurements  of  the  output 
of  each  of  the  transmitting  horns  were  made  in  the  anechoic  chamber 
while  the  carrier  and  modulation  frequencies  were  varied.  Figures  15 
to  18  are  examples  of  far-field  signals  measured  by  the  receiving  horn 
with  the  system  set  up  with  the  parameters  summarized  in  table  2.  The 
measured  far-field  time-  and  frequency-domain  output  from  the  carrier 
channel  (transmit  horn  1  in  fig.  14)  is  shown  in  figures  15  and  16;  the 
corresponding  time-  and  frequency-domain  output  from  the  modulation 
channel  (transmit  horn  2  in  fig.  14)  is  shown  in  figures  17  and  18.  The 
spectral  plots  clearly  show  that  the  IM  system  can  generate  pure  tones 
of  frequencies  fc  and  /c  —  fm  and  propagate  them  without  distortion 
into  free  space. 

To  verify  that  the  outputs  of  the  two  channels  of  the  IM  were 
matched  over  the  full  modulation  frequency  band,  we  swept  the  mod- 
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Table  2:  Sample  setup  parameters  for  far-field  IM  measurements. 


Parameter 

Value 

Carrier  frequency,  fc 

990  MHz 

Modulation  frequency,  fm 

463  kHz 

IM  sideband  suppressed 

Upper  sideband 

Transmit  horn  spacing 

43.2  cm 

Receive  horn  range 

430  cm 

Figure  15:  Measured  far-field  time-domain  signal  from  IM  carrier  chan¬ 
nel  (transmit  horn  1). 
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Figure  16:  Measured  far-field  frequency-domain  signal  from  IM  carrier 
channel  (transmit  horn  1). 


Figure  17:  Measured  far-field  time-domain  signal  from  IM  modulation 
channel  (transmit  horn  2). 
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fc  ~fm 


Frequency  (GHz) 

Figure  18;  Measured  far-field  frequency-domain  signal  from  IM  modu¬ 
lation  channel  (transmit  horn  2). 

ulation  input  from  100  Hz  to  1  MHz,  and  the  far-field  power  radiated 
from  each  channel  was  recorded  by  the  receiving  horn.  Figure  19  is  a 
plot  of  the  received  power  from  each  transmitting  horn  as  a  function 
of  the  modulating  frequency.  As  seen  in  the  figure,  the  amplitudes  of 
the  two  output  channels  are  well  matched  over  the  desired  frequency 
range. 

After  each  individual  channel  of  the  IM  system  was  evaluated,  the 
two  transmitting  horns  were  driven  simultaneously  (as  in  normal  oper¬ 
ation).  With  the  carrier  frequency  fixed  at  fc  =  990  MHz,  the  mod¬ 
ulating  frequency  was  swept  from  100  Hz  to  1  MHz  and  the  far-field 
time-  and  frequency-domain  behavior  of  the  IM  system  was  recorded 
by  the  receiving  antenna  (see  fig.  14).  Figures  20  to  23  are  plots  of  the 
measured  time-  and  frequency-domain  data  for  modulation  frequencies 
of  463  kHz  and  100  Hz.  Note  that  for  clarity,  the  abscissa  of  figure  23 
plots  the  measured  frequency,  /,  minus  the  carrier  frequency,  fc  =  990 
MHz.  As  can  be  seen  from  the  figures,  the  radiated  waveforms  from  the 
two  channels  of  the  IM  system  linearly  combine  in  the  far  field  to  pro¬ 
duce  a  waveform  with  a  duty  factor  of  50  percent,  a  carrier  frequency, 
fc,  equal  to  the  input  carrier  frequency,  and  a  modulation  envelope  fre¬ 
quency  equal  to  the  difference  between  the  input  carrier  and  modulation 
frequencies  {fc  -  fm)- 
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Figure  19:  Received  power  as  a  function  of  modulating  frequency,  f^, 
from  each  of  two  transmitting  horns. 


Figure  20:  Measured  far-field  (r  =  430  cm)  time-domain  behavior  of 
IM  for  /e  =  990  MHz,  U  =  463  kHz. 
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Figure  21:  Measured  far-field  (r  =  430  cm)  frequency-domain  behavior 
of  IM  for  fc  =  990  MHz,  fm  =  463  kHz. 


Figure  22:  Measured  far-field  (r  =  430  cm)  time-domain  behavior  of 
IM  for  fc  =  990  MHz,  =  100  Hz. 


Figure  23:  Measured  far-field  (r  =  430  cm)  frequency-domain  behavior 
of  IM  for  /c  =  990  MHz,  fm  =  100  Hz. 


4  Conclusions 


The  IM  is  a  broadband  microwave  source  with  independently  timable 
carrier  and  modulation  frequencies.  A  general  methodology  for  the 
design  of  an  n-channel  (n  >  2)  IM  has  been  developed,  along  with 
computer  codes  (see  the  appendices)  to  assist  in  the  design  process.  To 
validate  the  design  methodology,  we  developed  a  prototype  two-channel 
device.  We  give  details  of  the  design  and  fabrication  of  key  subsystems, 
such  as  a  broadband  90°  phase-difference  network  and  single-sideband 
modulator.  The  prototype  IM  was  experimentally  shown  to  have  met 
its  design  specifications  for  the  carrier  frequency  (0.8  <  /c  <  2  GHz) 
and  modulation  frequency  (100  Hz  <  /„,  <  1). 


28 


A  Code  Listing  for  Pole-Zero  Solver 


//  phaseshift.cc 

#def  ine  MAIN 

#include  ’’dblvec.h” 

#include  ’’dblply.h” 

#include  ”  cmplxdbl .  h” 
tinclude  ” graphic2d .  h” 

double  ac[5]  ={1.38629436112,0.09666344259,0.03590092383, 

0.03742563713,0.01451196212}; 

double  bc[5]  = 

{0.5,0.12498593597,0.06880248576,0.3328355346,0.00441787012}; 

DblPoly  Ka(4,ac); 

DblPoly  Kb(4,bc); 

double  Kelliptic(double  m); 

double  EllipticJlatio(double  q,double  r, double  tol  =  l.e-6,int 
nmax  =  20); 

CmplxDbl  Filter Jlesponse(double  w,DblVec  poles); 

//  Complete  Elliptic  Integral  of  the  1st  kind 

double  Kelliptic(double  m) 

{ 

double  ml  =  1.0-m; 
return(Ka(ml)+log(1.0/ml)*Kb(ml)); 

} 

//  Elliptic  Ratio 

double  Elliptic-Ratio (double  q,double  r,double  tol,int  nmax) 

{ 

double  num0,den0,rat0; 
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double  num  =  cos(r); 
double  den  =  sin(r); 
double  rat  =  num/den; 


double  qi  =  1.0; 
double  q2  =  q*q; 
double  qrat  =  1.0; 

for  (  int  i  =  1;  i  <  nmax;  ++i) 

{ 

qrat  *  =  q2; 
qi  *  =  qrat; 

ratO  =  rat; 
numO  =  num; 
denO  =  den; 

num  +=  qi*cos(((double)(2=t:i+l))*r); 

if  ( i%2  ^  0  ) 

{ 

den  -=  qi*sin(((double)(2*i+l))*r); 

} 

else 

{ 

den  +=  qi*sin(((double)(2*i+l))*r); 

} 

rat  =  num/den; 

if  (  fabs(rat-ratO)  <  tol*fabs(rat)  )  retinrn(rat); 

} 

error("\n7,d  iterations  not  enough  for  7,g 
t  olerance\n "  ,nmax,tol) ; 

return(O.O); 
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} 


//  Calculate  response  of  an  all  pass  filter  with  a  given  set  of  poles 

CmplxDbl  Filter-Response  (double  w,DblVec  poles) 

{ 

CmplxDbl  Z(1.0,0.0); 

CmplxDbl  z(0.0,w); 

for  ( int  i  =  0;  i  <  poles.dim();  ++i) 

{ 

Z  *  =  (z-poles(i))/(z+poles(i)); 

} 

return(Z); 

} 

main(int  argc,char  **argv) 

{ 

FILE  =t=outfile  =  NULL; 

if  (  argv[2]  NULL  )  outfile  =  createascii(argv[2]); 

ArgLst  *args  =  new  ArgLst(argv[l]); 


A 


Inputs  for  90  degree  phase  shift  filter  synthesis  program 


npoles  =  number  of  poles  in  filter  (default  8) 
flo  =  Minimum  frequency  (Hz)  at  which  90  degree 
phase  shift  is  required  (default  100) 
fhi  =  Maximum  frequency  (Hz)  at  which  90  degree 
phase  shift  is  required  (default  led) 
nplot  =  Number  of  frequency  points  to  calculate  filter 
response  at  (default  101) 

flo-plot  =  Minimum  plot  frequency  (Hz)  (default  10) 
fliLplot  =  Maximum  plot  frequency  (Hz)  (default  lOeS) 
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♦/ 


Int  *npoles_  =  new  Int("npoles",NULL,"8"); 
axgs^addarg((Arg  *)npoles_); 

Double  *flo_  =  new  Double("flo", NULL, "100"); 
args— >axldaxg((Arg  *)flo_); 

Double  ♦fhi_  =  new  Double("fhi", NULL, "1000000"); 
axgs— »^addarg((Arg  *)fhi_); 

Int  *nplot_  =  new  Int("nplot",NULL,"101"); 
args— »addarg((Arg  *)nplot_); 

Float  ♦flo_plot_  =  new  Float("flo_plot",NULL,"10"); 
args^addarg((Arg  *)flo_plot_); 

Float  *fhi_plot_  =  new  Float("fhi_plot",NULL,"  10000000"); 
args— »^addarg((Arg  *)fhi_plot_); 

args-^parse(); 

double  flo  =  flo_— >val(); 
double  fhi  =  fhi_— >val(); 
int  npoles  =  npoles_— >val(); 

//  Calculate  pole  and  zero  locations  for  filter  parameters 

double  Wrat  =  fhi/flo; 

double  WratJnv  =  1.0/ Wrat; 
double  Wrat_root  =  sqrt(Wrat); 
double  m  =  1.0- WratJnv* WratJnv; 
double  k  =  sqrt(m); 
double  kp  =  sqrt(l.O-m); 
double  theta  =  asin(k); 
double  ml  =  1.0-m; 
double  K  =  Kelliptic(m); 
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double  Kp  =  Kelliptic(ml); 

double  q  =  exp(-pi*Kp/K); 

double  delta  =  4.0=t=pow(q,(double)npoles); 

double  dB -Reject  =  20*logl0(1.0/tan(delta/2)); 

//  Print  locations  of  poles  and  zeros  for  circuit  synthesis 

printf("\n(fhi/flo)  =  7,lg  n  =  %d  flo  =  */,lg  Hz  fhi  =  7,lg 
Hz "  ,Wrat,npoles,flo,fhi); 

printf("\ndelta  =  %lg  deg  q  =  */i0.81g  theta  = 

70 . 81g\n",delta/dgtord,q, 
theta/dgtord); 

printf("\nSideband  Rejection  =  7«lg  dB\n'',dB -Reject); 

DblVec  pole(npoles); 

for  (  int  ipole  =  0;  ipole  <  npoles;  ++ipole) 

{ 

double  r  Jpole  = 

( (double)  ( 2  *ipole+ 1 ) )  *pi/  (double)  (4=i'npoles) ; 

if  ( ipole%2  ^  0  )  rJpole  =  -rJpole; 

pole(ipole)  =  Elliptic-Ratio(q, rJpole); 

printf("\nr7od  =  7.0.61g\t\t  p7.d  = 

70 . 61g", ipole, rJpole/dgtord, ipole, 
pole(ipole)); 

} 

cr(); 

float  flo-plot  =  flo-plot-— »^val(); 
float  fhi-plot  =  fhi_plot-^val(); 
int  nplot  =  nplot-— >val(); 

float  logJac  =  log(fhi-plot/flo-plot)/(float)nplot; 
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FltVec  £req(nplot); 

FltVec  phase(nplot); 

FltVec  amp(nplot); 

float  fiaorm  =  sqrt(flo*fhi); 

for  ( int  iplot  =  0;  iplot  <  nplot;  ++iplot) 

{ 

freq(iplot)  =  flo_plot*exp((float)iplot*log_fac); 
double  omega  =  (double)  (freq(iplot)/fhorm); 

CmplxDbl  Z  =  Filter_Response(omega,pole); 

phase(iplot)  =  Z.phasedeg(); 
amp(iplot)  =  Z.abs(); 

} 

//printcols(”  Hz  deg”, 

NULLjStdout,  1 024,  &freq,  &phase,  &aTnp,NULL ); 

//  Plot  response  of  filter 

Graph  *grph0  =  GraphJnit("Xll", "postscript", outfile); 
Graph  ♦grph.plt  =  new  Graph(grphO, 0.65, 0.75, 0.1, 0.1); 
Graph  *grph_title  =  new  Graph(grph0,0.65,0.1,0.1,0.9); 
Graph  *grph_poles  =  new  Graph(grph0,0.2,0.75,0.77,0.1); 

Axis  =t!Xax  =  new  Axis(0,0,0,"Hz"); 

Axis  *Yax  =  new  Axis(0,0,0, "degrees"); 

Xax— >setlog(); 

XYplot  *Phase_plot  =  new  XYplot(grph_plt,Xax,Yax); 

Phase_plot— >setframe() ; 

Phase-plot  ^setgrid(RED,NULL) ; 
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Curve2d  ♦phase.crv  =  new  Curve2d(&;freq,&phase); 

Text  Win  ♦title  =  new  Text  Win  (grph.title); 

title— >addtext(" Constant  90  deg  7,d  pole  Phase  Shift 
Filter  ",npoles); 

title^addtext(''flo  =  7,lg  Hz  fhi  =  7olg  Hz",flo,flii); 
title— >addtext("delta  =  7.1g  deg  Rejection  =  7.1g 
dB"  ,delta/dgtord,dB -Reject); 

Text  Win  ♦poles  =  new  TextWin(grph_poles); 

poles-^setxalign(LEFT -ALIGN) ; 
poles^addtext  ("Zeros"); 

double  freq-scal  =  sqrt(flo^fhi); 

for  ( ipole  =  0;  ipole  <  pole.dim();  ++ipole) 

{ 

poles— ♦addtext("  7+0 . 61e",pole(ipole)); 

} 

Phase.plot— >^addtext(title) ; 

Phase-plot— >^addtext(poles) ; 

Phase-plot-^addcurve(phase-crv) ; 

Phase-plot— »draw() ; 

Phase-plot^draw(PROMPT) ; 

Graph-term(); 

for  ( ipole  =  0;  ipole  <  pole.dim();  ipole  +=  2  ) 

{ 

double  p-plus  =  twopi^freq-ScaRpole(ipole); 
double  p_minus  =  twopi^fi:eq_scaRpole(ipole+2); 
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double  p-avg  =  (p_plus+p-ininus)/2.0; 
double  p_del  =  p_plus-p_miniis; 

printf("\np+  =  7.+061e  p-  =  7,+0 . 61e\n<p>  =  y.+O.Sle 
dp  =  7.+0.61e", 

p_plus,p_tninus,p_avg,p_del) ; 

double  D  =  sqrt((p_del*p_del-p_avg*p_avg)/3.0); 

double  tau-pliis  =  p-avg-D; 
double  tau_minus  =  (p_avg+D)/2.0; 

tau_plus  =  1.0/tau_plus; 
tau-ininus  =  1.0/tau_minus; 

double  K  =  2*tau_niinus/tau_plus+3; 

printf("\ntau+  =  7.+0.61e  tau-  =  7.+0.61e  K  =  = 
7.+0.61e", 

tau_plus,tau_minus,K) ; 


} 

cr(); 

exit(O); 


} 
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B  Code  Listing  for  Component  Value  Generator 

#include  <stdio.h> 

#include  <stdlib.h> 

#include  <math.h> 


#def ine  PI  3.1415926 
#define  TRUE  1 
#def  ine  FALSE  0 
#define  DimArrayl  7 
#define  MaxPoles  12 


/^H:**************************************************/ 


void  pause  0 

{ 

int  c; 

printf("Press  <return>  to  continue:  "); 
c=getchar(); 

} 


double  discrim  (pa, pb) 
double  pa,  pb; 

{ 

return(  (pa*pa)-(6*pa*pb) + (pb*pb) ) ; 

} 


void  synthesize_filter(n,wc,T) 
int  n; 

double  wc; 

double  T[DimArrayl]; 

{ 
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int  i,j,k,l,m,npair,iter; 
int  neg-flag, quit, swap; 
ch2U'  answer; 

double  temp,Q[DiraArrayl]; 
double  pa,pb; 

double  rlnrm,r2nrm,K,Xc,C,teinp2; 


/*  iterate  through  the  various  combinations  of  pole  pairs 
/*  load  the  pole  pair  set  into  the  array  QfiJ 


iter=l; 


for(i=l;i<(n-l);i++) 

{ 

Q[1]=T[1];  /*  first  pole  is  always  the  same  as 

Q[2]=T[i+l]; 

for(j=3;j<n;j++) 

{ 

Q[i]=T[j]; 

else 

QD]=T[2]; 

} 

for(k=l;k<(n-3);k++) 

/*  shuffle  the  rest  of  the  pole  pairs  */ 

{ 

switch  (k) 

{ 

case  1: 

{ 

breeik; 

} 

case  2: 

{ 

temp=Q[4]; 

Q[4]=Q[5]; 


*/ 

V 


Til]  */ 
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Q[5]=temp; 

bre2ik; 

} 

case  3: 

{ 

temp=Q[4]; 

Q[4]=Q[6]; 

Q[6]=temp; 

break; 

} 

default: 

{ 

printf("HEY!  I  got  nothin’  to  do!!!!\n"); 
break; 

} 

} 

printf( " \nlt erat ion  #%d :  \n"  ,iter) ; 
iter=iter+l; 

/*  print  out  the  particular  pole  pairs  */ 

for(l=l;l<(n-l);l++) 
printf("  y.g",  Q[l]); 
printf("  '/.gXn",  Q[n]); 

/*  test  to  see  if  any  discriminants  are  negative  */ 

negJiag=FALSE; 

for(m=l;m<(n-l);m=m+2) 

{ 

pa=Q[m]; 

pb=Q[in+l]; 

/*  printf(”  discrim[%d]  =  %g\n”  ,m,discrim(pa,pb));*/ 

if  (discrim(pa,pb)  <0)  neg_flag=TRUE; 


/*  if  no  discriminants  are  <  0,  then  synthesize  filter  */ 
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if(negJag^TRUE) 

{ 

printf("\n  Passed  discriminant  test...\n"); 

qiiit=FALSE; 

swap=FALSE; 

while(quit#TRUE) 

{ 

npair=l; 

for(m=l;m<(n-l);m=m+2) 

{ 

/*  compute  normalizing  reactance,  Xc  */ 

printf("\n  Input  C  (uF)  for  pole  pair  #y.d: 

",npair); 

scanf("%lf  ",  &:C); 

Xc=1.0/ (C*1.0e-6*wc); 
printfC  Xc  =  %g\n",  Xc); 

/*  compute  normalized  Rl  and  R2  */ 

pa=Q[ni]; 

pb=Q[in+l]; 

r2nrm=fabs(((pa+pb)  + 
sqrt(discrim(pa,pb)))/(2*pa*pb)); 

rlnrm=fabs(1.0/  (pa*pb=('r2nrm)); 

/*  if  swap=TRUE,  do  Rl  <-  R2/2  and  R2  <- 
2R1  to  try  to  get  */ 

/*  more  favorable  K  value 

if(swap==TRUE) 

{ 

temp2=r2nrm; 

r2nrm=2*rlnrm; 
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rlnrm=temp2/2; 


} 

K=3+ (2*r2iirm /rlnrm) ; 


/k.  pnntf(”\n  Rlnrm  =  %g 

R2nrm  =  %g\n”,rlnrm,r2nrm);*/ 

printf(’'\n  R1  =  7.g  R2  =  ‘/.g  K  =  ’/.gV", 
(r  lnrm*Xc) ,  (r2nrin*Xc)  ,K) ; 
npair=npair+l; 

} 


printf("\n  Try  again  <y/n>?  "); 

scanf(  "*/,s  "  ,&answer) ; 
if  ((answer==’y’)  [  |  (answer==’Y’)) 

{ 

swap=FALSE; 

printf("\n  Swap  R1  and  R2  <y/n>?  "); 

scanf("%s",  &answer); 
if  ( (answer==’y ’)  1 1  (answer==’ Y’) ) 
swap=TRUE; 

} 

else  quit=TRUE; 

} 

} 

else  printf("\n  Failed  discriminant  test  (discrim  < 

0)...\n"); 

} 

} 

} 


main() 

{ 

FILE  *infile; 

int  i,j,lim,npoles; 

double  temp,P[MaxPoles+l],S[7]; 
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double  fa,  fb,  wa,  wb,  wc; 


/*  get  total  number  of  poles  and  pole  values  from  input  file  */ 

infile=fopen( " input . dat ",  " r " ) ; 
fscanf(infile,  &;npoles); 

for(i=l;i<npoles;i++) 

fscanf(infile,  "7,lf ",  &;P[i]); 

/*  get  frequency  limits  */ 

fscanf(infile,  "7,lf/,lf ",  &fa,  &:fb); 

wa=2*PI*fa; 

wb=2*PI*fb; 

wc=sqrt  ( wa*wb) ; 

fclose  (infile); 

/*  order  poles  from  most  positive  to  least  positive  */ 

/*  and  store  in  P[ij  ♦/ 

i=l; 

while(i<npoles) 

{ 

for(j=(i+l);j<npoles;j++) 

{ 

if(P[i]  <  P[j]) 

{ 

temp=P[i]; 

P[i]=P[j]; 

P[j]=temp; 

} 

} 

i=i+l; 

} 


/*  load  in  positive  valued  poles  */ 
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for(i=l;i<(npoles/2);i++) 

S[i]=P[i]; 

printf("\n  PROCESSING  POSITIVE  VALUED  ROOTS  (I)\n"); 
printf(" 

synthesize-filter  ((npoles  /  2)  ,wc,S) ; 


/*  load  in  negative  valued  poles  */ 


for(i=l;i<(npoles/2);i++) 

S[i]=P[i+(npoles/2)]; 


printf("\n  PROCESSING  NEGATIVE  VALUED  ROOTS  (Q)\n"); 
printf(" 


synthesizeJilter(  (npoles /2) ,  wc,S) ; 


} 
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C  SPICE  Input  Files 


C.l  90°  Phase-Difference  Network  —  Ideal  Op-Amp  Model 

3|c  *  *  ak  *  9|e  9(c  a|e  ^  3(e  *  9tc  9k  9|e  ^ 

♦ 

*  12  pole  I  &  Q  all-pass  network 

*  fa=100  Hz  fb=l  MHz  DKA  3/8/94 

*  ideal  op  amp  model 

* 

9|e  :k  ak  9k  9f(  afc  ak  ak  a|e  ak  ak  9k  9|e  9k  3k  ak  ak  ak  ♦  ak  9k  9|c  ak  9k  3k  9k  ak  9k  9k  9k 

vin  1  0  0.0  ac  1.0  0.0 

3k  3k  9k  9k  9k  ak  ak  9k  ak  ak  ak  ak  ak  9k  ak  9k  3k  ak  3k  9k  ak  ak  ak  ak  3k  3k  9k  ak  ak  9k  3k 

*  op  amp  1  (I  ckt) 

3k  ak  9k  3k  ak  ak  ak  ak  3k  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  9k  ak  ak  ak  9k  ak  9k 

xampl  325  opamp 
rll  1  3  310.0 

cll  1  3  O.OOle-6 

rl2  3  4  167.0 

cl2  4  0  O.OOle-6 

rl3  1  2  1000.0 

rl4  25  4079.0 

3k  ak  3k  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  9k  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  9k  ak  ak 

*  op  amp  2  (I  ckt) 

3k  3k  9k  3k  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  9k  9k 

xamp2  769  opamp 
r21  5  7  762.0 

c21  5  7  0. Ole-6 

r22  7  8  1374.0 

c22  8  0  0. Ole-6 

r23  5  6  1000.0 

r24  6  9  6608.0 

9k  9k  9k  ak  9k  3k  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak  ak 

*  op  amp  3  (I  ckt) 

:|c  :|c  *  3|c  :|c :|c:tc  4:  %  :|c  *  :|c  :|c  itc  9|c  :tc  :te  *  !(e  4;  9|(  :|e  3f;  itt  4c 

xamp3  11  10  13  opamp 

r31  9  11  3242.0 

c31  9  11  O.le-6 

r32  11  12  4075.0 
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c32 

12 

0  O.le-6 

r33 

9 

10  1000.0 

r34 

10 

13  5514.0 

rl3 

13 

0  50.0 

*  5|c  :|e  *  :|c  ♦  3|e  :^c  *♦*♦*♦*♦♦  3|c  3|e  3|e  ♦  :4c  s|c  :*c  3|e  j|c  3»c  * 

*  op 

amp 

4  (Q  ckt) 

*  :4c  *  :|c  ajc  *  *  ♦  :)c  *  :4c  jjc  *  *  :|c  *  ♦  ♦  :4c  3f:  :4c  :►  :4c  *  :4c  ♦  afc  ♦  ♦  :4e 

xamp4  23 

22  25  opamp 

r41 

1 

23  7569.0 

c41 

1 

23  0.2e-6 

r42 

23 

24  4084.0 

c42 

24 

0  0.2e-6 

r43 

1 

22  1000.0 

r44 

22 

25  4079.0 

**:4c*:4c*5|c************************ 

♦  op 

amp 

5  (Q  ckt) 

:4c*****:|c***^*ic*%:4c****:4c*34s****3|cs4c*34c 

xamp5  27 

26  29  opamp 

r51 

25 

27  1843.0 

c51 

25 

27  0. Ole-6 

r52 

27 

28  3325.0 

c52 

28 

0  0. Ole-6 

r53 

25 

26  1000.0 

r54 

26 

29  6608.0 

:4c  :4c  *  ♦  *  *  ♦  :4c  :4c  :4c  afc  *  :4c  *  ♦  ♦  :4c  :4c  :4c  :4c  :4c  *  *  *  ♦  *  :4c  :4c  ♦  :4c  :4c 

*  op 

amp 

6  (Q  ckt) 

:4c  :4c  :4c  :f:  :4c  :4c  :4c  :4c  :4;  :4c  :4c  :4c  :4c  :4c  :4c  :4c  :4c  :4c  :4c  :4c  :4c  :4c  :4c  :4c  :4c  %  :4c  :4c  :4c  :4c  :4c 

xamp6  31 

30  33  opamp 

r61 

29 

31  622.0 

c61 

29 

31  O.OOle-6 

r62 

31 

32  781.0 

c62 

32 

0  O.OOle-6 

r63 

29 

30  1000.0 

r64 

30 

33  5514.0 

rl6 

33 

0  50.0 

af:  3|c  *  j|e  ♦  ajcjle**  a(c  **♦*♦♦*  s|e  5|s  :►*  +  **♦♦*  j|e 

*  ideal  op  amp  w/  100k  gain  and 

*  one-pole  roll  off  at  approx.  1  MHz 
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:|c  4c  :|c  :(e  afe  :(e  9^  :|e  :(e  9|e  :tc  3tc ;tc:|c  3|c  :|c  :|e  4:  %  aie  *  3|c 

.subckt  opamp  non  inv  out 

rin  non  inv  1.0e6 

egain  1  0  non  inv  100. 0e3 

ropen  1  2  5.0e3 

copen  2  0  30.0e-12 

eout  30  201 

rout  3  out  50 

.ends  opamp 

4c  a|e  ate  ak  ♦  a|e  *  4^  4e  4c  4e  4c  4c  afc  4e  *  3|c  4c  4e  4e  4c  4c  *  3ic  4c  4e  4e  4e  4c  3|c 

.ac  dec  100  1.0  3e6 
.end 
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C.2  90°  Phase-Difference  Network  —  AD846  OpnAmp  Model 

:|e  *  3(c  9k  :|c  :|e  3|o|e  3|e  ie  3^  :4c 

♦ 

*  12  pole  I  &  Q  all-pass  network 

*  fa=100  Hz  fb=l  MHz  DKA  4/19/94 

*  Analog  Devices  AD846  op  amp  model; 

*  "real"  chip  R  ft  C  values 

*  (same  basic  circuit  as  IQ2C_12P.CIR 

*  with  some  diagnostic  additions  and 

*  output  attenuator  networks) 

* 

*  I  and  Q  outputs  drive  50  ohms 

* 

*  Q  network  has  0.95  attenuation  on  output 

♦ 

9k***5|e**:k**  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

vin  1  0  0.0  ac  1.4  0.0 
vpos  100  0  dc  18.0 
vneg  101  0  dc  -18.0 

*  vsl  and  vs2  are  pos.  and  neg. 

*  supply  current  monitors 
vsl  100  90 

vs2  101  91 

*  3tc :tc  *  :tc  :tc  *  3):  lie**  ^  iK  ***:)(  :te  i|c  ic  ic  *  3|c  3)c 

*  op  amp  1  (I  ckt) 

******************************* 
xampl  3  2  90  91  5  50  AD846 

rll  1  3  309.0 

cll  1  3  O.OOle-6 

rl2  3  4  169.0 

cl2  4  0  O.OOle-6 

rl3  1  2  1000.0 

rl4  2  5  4120.0 

9k:tc:(c9k:k3k:(c3|c:(c9k:k3fc:|c9k:k:tc:|c9k:k:|c:|c:fc:4c9k4c:(c9k:}c9ka|e:|c 

*  op  amp  2  (I  ckt) 

3k  :k  3k  :fc  3k  :k  :4c  :k  :|c  9tc  :4c  9k  :|c  9k  :k  :|c  9k  :|C  3k  :4c  3k  :|c  :fc  9k  :k  3k  3k  :|c  9|C  :fc  9k 

xamp2  7  6  90  91  9  51  AD846 
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r21 

5 

7 

768.0 

c21 

5 

7 

0.  Ole-6 

r22 

7 

8 

1370.0 

c22 

8 

0 

0.  Ole-6 

r23 

5 

6 

1000.0 

r24 

6 

9 

6650.0 

afe  af;  :4c  3|e  a|c  a|c  a|e  a|c  a|c  a(e  3^  ale  a|e  a|e  :|e  a^  afe  3|c  a(e  a|c  a(e  ate  a(e  af(  a|e  a|e  a|c  ate  a|c  a|e  a)e 

*  op 

amp 

3 

(I  ckt) 

afe  aje  ate  a4e  3|e  a|e  ate  a|e  a|e  a|e  a|e  a|c  a(e  ate  a|c  a|e  ate  a)e  afe  a|e  a»;  a)e  a|e  a|e  af;  a|e  ate  ate  ate  ate  a|e 

xamp3  11 

10 

90  91  13  52  AD846 

r31 

9 

11 

3240.0 

c31 

9 

11 

O.le-6 

r32 

11 

12 

4120.0 

c32 

12 

0 

O.le-6 

r33 

9 

10 

1000.0 

r34 

10 

13 

5490.0 

rs3 

13 

14 

105.0 

rp3 

14 

0 

2.0e3 

rl3 

14 

0 

50.0 

a|e  ate  ate  ate  ate  ate  ate  a|e  ate  ate  ate  ate  ate  ate  a|e  ate  ate  ate  ate  ate  ate  ate  a|e  ate  ate  ate  a|e  a|e  ate  ate  ate 

♦  op 

amp 

4  1 

CQ  ckt) 

ate  afe  ate  ate  a|e  ate  ate  ate  ate  ate  ate  ate  ate  ate  ate  ate  ate  ate  afe  ate  a|e  ate  ate  ate  ate  ate  ate  ate  ate  ate  ate 

xamp4  23  22 

90  91  25  53  AD846 

r41 

1 

23 

7500.0 

c41 

1 

23 

0 . 2e— 6 

r42 

23 

24 

4120.0 

c42 

24 

0 

0.2e-6 

r43 

1 

22 

1000.0 

r44 

22 

25 

4120.0 

*  op  amp  5  (Q  ckt) 

:te  *  *  :|c :(( 3fc  9|e  :|c  9|c :(( :fc  3|e  afc  3(e  9|c  :|c  3)c  aK  Ne  9fe  *  3|c  :|e  ate  3fc  *  3fe  ate  :te  :te 

xamp5  27  26  90  91  29  54  AD846 

r51  25  27  1820.0 

c51  25  27  0. Ole-6 

r52  27  28  3320.0 

c52  28  0  0. Ole-6 

r53  25  26  1000.0 
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r54  26  29  6650.0 

:|c  4;  :)c  I):  4:  *  *  4=  *  *  *  *  1)1  *  *  *  :•=  *  =4^  4c  4:  *  *  *  *  :^  *  *  *  * 

*  op  amp  6  (Q  ckt) 


xamp6  31 

30 

90  91  33  55  AD846 

r61 

29 

31 

619.0 

c61 

29 

31 

O.OOle- 

-6 

r62 

31 

32 

787.0 

c62 

32 

0 

O.OOle- 

-6 

r63 

29 

30 

1000.0 

r64 

30 

33 

5490.0 

rs6 

33 

34 

105.0 

rp6 

34 

0 

2.0e3 

rl6 

34 

0 

50.0 

3»c  j|e  a|c  *  :Jc  :|c  a|c  :►  3»£  ****:►♦♦♦**♦♦♦  :|c  3f£  :fc  + 
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*  Node  assignments 


* 

* 

* 

* 

* 

* 

* 

.SUBCKT  AD846 
♦ 

♦  INPUT  STAGE 

♦ 

R1  99  8 


non- inverting  input 
inverting  input 
positive  supply 
I  negative  supply 
I  I  output 

I  I  I  compensation  pin 

till 

1  2  99  50  28  12 


1E3 
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R2 

10 

50 

1E3 

VI 

99 

9 

1.7 

D1 

9 

8 

DX 

V2 

11 

50 

1.7 

D2 

10 

11 

DX 

11 

99 

5 

275E-6 

12 

4 

50 

275E-6 

Ql 

50 

3 

5  QP 

Q2 

99 

3 

4  QN 

Q3 

8 

6 

30  QN 

Q4 

10 

7 

30  QP 

R3 

5 

6 

300E3  - 

R4 

4 

7 

300E3 

* 

*  INPUT 

ERROR  SOURCES 

* 

BGl 

99 

1 

I=3E-6  +  (5E-9  * 

V(l,22)) 

BG2 

99 

30 

I=150E-9  +  (5E-9 

♦  V(l,22)) 

VOS 

3 

1 

25E-6 

LSI 

30 

2 

50E-9 

CSl 

99 

2 

2.0E-12 

CS2 

99 

1 

2.0E-12 

♦ 

EREF 

97 

0 

22  0  1 

* 

♦  TRANSCONDUCTANCE  STAGE 

♦ 

R5 

12 

97 

200E6 

C3 

12 

97 

5.5E-12 

G1 

97 

12 

99  8  lE-3 

G2 

12 

97 

10  50  lE-3 

V3 

99 

13 

1.3 

V4 

14 

50 

1.3 

D3 

12 

13 

DX 

D4 

14 

12 

DX 

♦ 


♦  POLE  AT  100  MHZ 

♦ 


50 


R8 

17 

97 

1E6 

C4 

17 

97 

1.59E-15 

G4 

97 

17 

12  22  lE-6 

♦ 

♦  POLE 

AT 

100  MHZ 

♦ 

R9 

18 

97 

1E6 

C5 

18 

97 

1.59E-15 

G5 

97 

18 

17  22  lE-6 

* 

*  POLE  AT  200  MHZ 

* 

RIO  19  97  1E6 

C6  19  97  0.795E-15 

G6  97  19  18  22  lE-6 

♦ 

*  POLE  AT  500  MHZ 

♦ 

Rll  20  97  1E6 

C7  20  97  0.318E-15 

G7  97  20  19  22  lE-6 

♦ 

♦  POLE  AT  500  MHZ 

* 

R12  21  97  1E6 

C8  21  97  0.318E-15 

G8  97  21  20  22  lE-6 

♦ 

*  OUTPUT  STAGE 

* 

ISY  99  50  3.55E-3 

R13  22  99  16.67E3 

R14  22  50  16.67E3 

R15  27  99  60 

R16  27  50  60 

L2  27  28  lE-8 

G9  25  50  21  27  16.67E-3 

GIO  26  50  27  21  16.67E-3 
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Gll  27  99  99  21  16.67E-3 

G12  50  27  21  50  16.67E-3 

V5  23  27  1.35 

V6  27  24  1.35 

D5  21  23  DX 

D6  24  21  DX 

D7  99  25  DX 

D8  99  26  DX 

D9  50  25  DY 

DIO  50  26  DY 

* 

♦  MODELS  USED 

* 

.MODEL  QN  NPN(BF=1E9  IS=1E-15  VAF=16.3E3) 

.MODEL  QP  PNP(BF*1E9  IS=1E-15  VAF=16.3E3) 

.MODEL  DX  D(IS=1E-15) 

.MODEL  DY  D(IS=1E-15  BV=50) 

.ENDS  AD846 

**************^cH:^:*^i*********m** 

.ac  dec  100  1.0  3e6 
.end 
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